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FOREWORD 

This research was supported under contract NAS1-1Q277 with the National 
Aeronautics and Space Administration's Langley Research Center. The technical 
monitor on the contract was H. Harris Hamilton XX. The results are published 
in two parts: 

Fart I Description of Basic Method (NASA CR 111921) 

Part XI Description of Computer Program (NASA CR 111922) 

Part I contains a detailed description of the theoretical approach and basic 
equations that were used to obtain a solution to the problem. Part II contains 
a detailed description o;f the computer program and is intended to serve 
primarily as a users' manual. 
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CALCULATION OF INVISCIL SURFACE 

streamlines and heat transfer 

on SHUTTLE TYRE CONFIGURATIONS 
I _ Description of Ba.lt Method 


Part 


By Fred R. DeJarnette 
North Carolina Star. University 
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INTRODUCTION 

Space shuttle vehicles ere subjected to appropriate thermal 

tb . reentry phase of their trajectory • J” „ p0B . d surfaces are - 

'protection systems. and flight conditions are being 

quired. Since many differen 


considered for the space shuttle, existing experimental heating rates are in- 
sufficient, and an exhaustive new experimental program would he very costly. 

This part of the report presents an analytical method for calculating laminar, 
transitional, and turbulent heating rates on general three-dimensional bodies 
for a wide range of reentry flight conditions. A description of the computer 
program appears in Part II (NASA CR 111922) . 

For those flight conditions where aerodynamic heating is Important, the 
flow field generally can be represented by a boundary layer adjacent to the 
surface and an inviscid layer between the boundary layer and bow shock wave. 
Reentry velocities are high enough to cause the gas in the shock layer to devi- 
ate from a perfect gas, but sufficiently low (less than 30,000 ft/sec) that 
radiative heat transfer is small compared to convective heating. 

Numerical methods are currently available for calculating inviscid and 
boundary-layer flow-fields over two-dimensional and axisymmetric bodies at zero 
angle-of-attack (e.g. refs. 1 and 2). Also, methods for calculating inviscid 
flow fields over pointed and spherically blunted three-dimensional bodies have 
been developed (e.g. refs. 3 and 4). However, these techniques require con- 
siderable storage and computational time on existing digital computers, and a 
satisfactory three-dimensional boundary- layer program is not available. The 
problem is further complicated by the fact that shuttle type configurations are 
general three-dimensional bodies whose geometry is difficult to describe ana- 
lytically. During reentry they operate at large angles of attack, and laminar, 
transitional, and turbulent heating may occur. 

If tractable solutions of three-dimensional flow fields are to be obtained, 
simplifying approximations must be used for the inviscid and viscous flow- 
fields. A substantial simplification to the viscous flow field may be achieved 
by the axisymmetric analogue, or small cross flow approximation, for three- 
dimensional boundary layers (refs. 5 and 6). This approximation allowB the 
heat transfer rate along a surface inviscid streamline to be calculated by any 
method applicable to an equivalent body of revolution at zero angle-of-attack 
provided the surface inviscid solution is known. Along each streamline the 
"equivalent radius" or scale factor replaces the radius of the body of revo- 
lution in the axisymmetric analogue, and the distance along a streamline re- 
places distance along the body of revolution. Thus each surface inviscid 
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streamline on the three-dimensional body corresponds to a different equivalent 
body of revolution at zero angle-of-attack. 

In order to apply the axisymmetric analogue, the invlp:ld solution at the 
edge of the boundary layer must be known- For blunt nosed bodies with thin 
boundary layers, the flow field properties at the edge of the boundary layer 
are nearly the same as the corresponding lnvlscld solution on the surface- If 
the surface streamlines are assumed to pass through the normal part of the bow 
shock wave, then the surface entropy (for a perfect or equilibrium gas) is 
equal to the entropy aft of the normal shock wave. The remaining inviscid 
solution needed for the boundary layer solution can be calculated from the 
pressure distribution; and when the pressure distribution is unknown, modified 
Newtonian pressures can be used. 

One of the major difficulties in applying the axisymmetric analogue is the 
calculation of the surface inviscid streamlines and the corresponding "equiva- 
lent radius" or scale factor. References 7, 8, and 9 present several approach- 
es for calculating the surface inviscid streamlines. Approximate methods for 
determining the "equivalent radius" are given by Timmer (ref. 9), Pinkus and 
Cousin (ref. 10) and Vaglio-Laurin (ref. 11), whereas references 12 and 13 
determined equivalent radii on blunted cones at an angle of attack by calcu- 
lating the spacing between streamlines. These approximate techniques are not; 
satisfactory for general three-dimensional bodies. 

DeJarnette and Davis (ref. 14) determined an analytical expression for the 
"equivalent radius" on blunted cones at an angle of attack based on a simpll- 
fled method for calculating the streamlines. Later DeJarnette and Tai (ref. 
15) developed an "exact" method for calculating the streamlines and "equivalent 
radius" from surface pressure distributions on axisymmetric bodies at an angle 
of attack. Both techniques are extended to general three-dimensional bodies) in 
this report. 


+ This assumption is invalid in those regions where the boundary layer has 
"swallowed" the entropy layer. 

* DeJarnette also obtained analytical results for elliptic cones at an angle 
of attack, using this simplified method, while participating in the NASA- 
ASEE Summer Faculty Research Program at Ames Research Center, NASA, Summer, 
1969 (unpublished). 
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Another difficulty associated with the surface streamlines is their be- 
havior in a three-dimensional stagnation region. The stagnation point is a 
singular point (nodal point) , and thus numerical problems are normally encoun- 
tered there, These problems are solved in the present method by developing 
analytical expressions for the streamline geometry and "equivalent radius" in 
the neighborhood of the stagnation point ■ 

In tho application of the axisymmetric analogue along each surface invis- 
cid streamline, methods are required for the calculation of laminar, transi- 
tional, and : turbulent heating rates on an equivalent body of revolution at zero 
angle-of-attack. Techniques which involve the numerical solution of the appro- 
priate boundary layer equations, a.g. ref. 2, could be used; however, this 
approach would destroy the simplicity and short computational time desired in 
the present- method. Therefore, approximate solutions of the boundary layer 
equations cn a body of revolution wi-1-1- be applied, 

Lees’ method (ref. 16) was used with the axisymmetric analogue- to calcu- 
late laminar heating rates on blunted cones at an angle of attack in references 
14 and 15, -and the results compared well with experimental data. However, 

Lees' method neglects pressure gradient effects and is valid only for cold 
walls,. As .a better approximation, the present method uses the local similarity 
method of Beckwith and Cohen (ref. 17), which includes the effects of pressure 
gradients and nori-cold walls on laminar beating rates. However, special con- 
sideration! must be given to the applicat ion of thiB method to three-dimensional 
stagnation points. 

The analytical calculation of transitional heating rates and location of 
the transition region has been unsuccessful to date, and the reader is referred 
"to ref 2 for a discussion of this phenomenon , For the present work, transi- 
tional beating rates are computed as a weighted sum of the laminar and turbu- 
lent heating rates, using the weighting distribution of Dhawan and Narasimha 
(ref. 16)-. Along a surface inviscid streamline, the beginning and end of 
transition may be specified optionally by one of the following three param- 
eters 5 ; 

1. geometric location, or 

2. value of local momentum thickness Reynolds number, or 

3. value of integrated unit Reynolds number along a streamline. 
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Many methods have been developed for calculating turbulent boundary layer 
properties (see discussion In ref. 2). As a compromise between accuracy and 
simplicity, the technique developed by Spalding and Chi (ref. 19) for calcu- 
lating turbulent skin friction coefficients based on momentum thickness 
Reynolds number is used here. The momentum thickness is calculated from the 
integral method of Reshotko and Tucker (ref. 20) for use in Spalding and Chi's 
equations. Then the skin friction coefficients are converted to turbulent 
heating rates through von Karman's form of Reynolds analogy (see ref. 21). 

For the applications herein, both a perfect gas and equilibrium air are 
considered. Therefore, modifications to the methods of Spalding and Chi (ref. 
19), and Reshotko and Tucker (ref. 20) are needed to make them applicable to 
equilibrium air. To keep the computations down to a minimum, the correlation 
formulas of reference 22 are used for equilibrium air properties. 

Results from the computer program are presented for a blunted 15° half- 
angle cone, a blunted 70° delta wing, and the HL-10 lifting body as example 
applications of the basic method. 
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SYMBOLS 

parameters defined by eqs. (150) and (151) 
speed of Bound at edge of boundary layer, ft/sec 
constants used in eq. (11) 

coefficients for spline function defined by eq. (25) 

Ratio of principal velocity gradients, see eq. (107) 
coefficients for spline function defined by eq. (25) 
parameter used in eq. (95) 

2 

local skin friction coefficient, 2 t /(p U ) 

W © 6 

incompressible skin friction coefficient 
derivative along a streamline 

unit vectors on body surface given by eqs. (33), (35), and (36) 
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unit vectors in streamline coordinate system, eqs. (37), (34), 
and (38) 

unit vectors in cylindrical coordinate system, see fig- 3 

eccentricity of conic section 

body radius, measured from longitudinal axis 

cubic spline function for body radius 

function defined by eq, 0?7) 

function defined by eq- (171) 

function defined by eq- (29) 

scale factor in 8-direction, dq * h d8 

2 2 

adiabatic wail enthalpy, ft / sec 

2 2 

enthalpy inside and at edge of boundary layer, respectively, ft /sec 
2 2 

wall enthalpy, ft /Bee 

8 2 2 

reference enthalpy, 2»119 x 10 ft / sec 
scale factor in s-direction, ds ■ h d£ 

.compressible form factor, H • 6*/8 m 
incompressible form factor 
transformed form factor 

2 2 

static enthalpy aft of normal shock, ft /sec 

2 2 

free -stream stagnation enthalpy, ft /eec 
coefficient used in cubic spline function 
Mach number 

straight-line distance normal to body surface 

compressible and incompressible Stanton numbers, respectively 
2 

pressure, lb/ft 

parameters defined by eq. (12) 
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x,y ,2 


Prondtl number 

distance normal to streamline on body surface, dq - h dg 

2 

heat-transfer rate at wall, BTU/ft -sec 

oody radius (r • f) 

body radius of curvature at nose 

principal body radii of curvature at stagnation point 
momentum thickness Reynolds number 


Integrated unit Reynolds number, 


s p U dS 
K e e 


w ® 2 2 

gas constant for undissociated air, 1716 ft /sec -R 
position vector for points on body 

distance along a streamline, measured from stagnation point 
(dS ■ h g dO 

body coordinates at stagnation point 

entropy 

time, sec 

enthalpy ratio, h e /H g 
temperature, ®R 

velocity insa.de and at edge of boundary layer, respectively, ft/sec 

function defined by eq. (174) 

velocity aft of normal shock wave, ft/sec 

inviscid speed on surface, ft/sec 

inviscid velocity on surface, ft/sec 

velocity components along and S_, ft/sec 

unit vector in direction of freestream velocity vector 

weighting function given by eq. (161) 


Cartesian coordinates 
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'w 

5 ! 

w 


m 


•tr 


w 


angle o£ attack 

coordinate noxmel to streamline on body «*“» 
pressure gradient parameter given by eg. (W) 
ratio of specific heats 
body angle defined by eq. (30) 
displacement thickness 

body angle defined by eq. (3D 

radius of circle around stagnation point, 

wall enthalpy ratio, 
enthalpy gradient normal to well 

inclination engie of surface invi.cid 

momentum thickness, ft 

petemeter defined by eg. <163) 

coefficient of viscosity, slug/ft-sec 

coordinate along a .tr.amlin. <dS - b.«> 

parameter defined by eg. U«> 

3 

mass density, slug/ ft 

angle defined by eq. (47) 

2 

wall shear stress, lb /ft 

circumferential angle, see fig. 4 

angle defined by eq. (78) 


Subscripts: 


lam 

B 


edge of boundary layer 
laminar value 
stagnation point 
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curb turbulent value 

w evaluated at wall 

e evaluated on e-circle, (see fig. 6) 

» freestream conditions 

tre end of transition region 

tri beginning of transition region 

2 conditions aft of normal shock wave 

ANALYSIS 


Geometry for General Three-Dimensional Bodies 

Before heating rates can be calculated, a description of the body geometry 
is needed which yields slopes and radii of curvature in both the axial and cir- 
cumferential directions. If the body could be described by analytic equations, 
these results would be easily obtainable. However, the geometry of shuttle- 
type configurations cannot, in general, be represented by an analytic equation. 
Often only the coordinates of the body at several axial .stations are known, and 
the body geometry must be generated from this data. Therefore, a numerical 
method is required to generate the body shape and calculate slopes and radii of 
curvature at any position on the body. This method must be capable of perform- 
ing numerical interpolation and differentiation reasonably accurately. 

It is well known that numerical differentiation is generally inaccurate. 
However, the method of splines (see ref. 23) has been proven to be an affective 
method for numerical differentiation, integration, and interpolation. The 
effectiveness of the cubic spline to perform these operations accurately is a 
consequence of the strong convergence properties it possesses. Before discuss- 
ing the two-dimensional cubic spline, which is needed to describe a three-di- 
mensional body, the one-dimensional spline will be developed first. 

To illustrate the cubic spline for one dimension (one Independent varia- 
ble) , consider the coordinates of a body cross-section (viewed from the rear) 
with one plane of symmetry as shown in figure 1. 
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Figure 1. - Coordinates of body cross section (viewed from rear) 


For the interval 0 £ 4> £ it, a mesh A of body positions is given by (y^ » 
with i ® 1, 2, .... i . It is convenient to transform these coordinates into 
polar coordinates by the transformation 


where 




+ , J 2] 


1/2 


« tan" 1 (z^/y^) 


0 » <t>, < d> 2 < ... ■ it 

J m 


(1) 

( 2 ) 


and (|> is the independent variable. Then the cubic spline for this mesh yields 
a function which is (1) continuous with its first and second derivatives 
on (0, it], (2) coincides with a cubic in each subinterval S ♦ S ♦j 0 “ 

2, 3, ...» j^), and (3) satisfies f^(<t>j) “ f j (3 “ *•» 2 > 3 m ) • Note that 
the points 0 O , ..., ♦. need not be equally spaced. 



d 2f A U ) 

Denoting f"(<f>.) ■ — s- 4 — by m , the cubic spline on the Qubinterval 

A j j 

[ (j*.] is assumed to have a linear variation of the second derivative, i.e 
<3 i j 

(•flj, “ (<t> ~ 


f A ( ^ “ m j-l A + "j 




<3) 


where A^ -4^-4^ l* Then inte 8 ratill S the above expression and evaluating 
the constants of integration, the first derivative and spline functions on 
*jl are 


( 4 >, - *) 

£’(4>) “-m. , — n~TZ * 

A j-1 2 A<p A j 
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( + “ ♦, x ) 
Tl$ ~ 


(£. - £ 


1 zL 


A<t>j 


- <» 3 - vi> 


A<j> 


( 4 ) 


and 








(<f>j “ 4>) 

(<t> “ 

, VlV 

(<^ 1 “ ♦) 

1 6 A^j 

' + m. , 7 “ t 

j o A*j 

_ 4-1 6 

^j 

A . 2_ 

U “ 
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A^ 




(5) 


In order to apply the equations above, it remains to calculate the second 

derivatives m^ at the mesh points. By virtue of equations (3) and (5), the 

functions f"(4>) and f.(4>) are continuous on [0, it]. By equating the one-sided 
A A « 4- 

limits of the derivatives at each interior mesh point, f ^ (4> j ) “ )» the 

function f ^ becomes continuous on [0, n]. These relations yield 


lit 


m 


j-i 


+ (A4»j + A$ 


p j+i 


> i. + fi+ili+i . E - 

' O T A A A 




A 4 > 


j+1 


A(|> 


(6) 


j 


for j - 2, 3, . . . , J m - 1. 


Equation (6) gives (j - 2) simultaneous equations in the unknowns 

in. itu . Therefore two additional conditions must be specified, the "end 

2 Jm 

conditions", to determine all the m^'s. If the slopes at the end points, 
f*(0 + ) and f^(ir~), are known, the one-sided derivatives from equation (4) 
provide the relations 


2m^ + m2 


1 

6 ! 

*2 “ f l 
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A$ x 
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X^* 1 ' 3m linear "equations 

„ This is e system of J„ 
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Then by successive elimination of ., M n. - 1 from the 2nd, 3rd 

IU 

j^-th equations of (12), the following equivalent equation system results: 


" q j m J+l + U j 

(J “ J “1» J -2, ...» 1) 

itt tn 
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These relations yield m. , m ,, ..., m. successively. Once the m. 's are 

J m ^ 


'm 


evaluated, the functions f^($), f ^ (d>) » and f ^(<J>) can be easily calculated from 


equations (5), (4), and (3), respectively, for any value of 4> in the interval 
0 £ $ < ir. 

The three-dimensional bodies considered herein will be described by r ■ 
f(x,$) where x is distance along the axlB from the nose, $ is the circumferen- 
tial position, (j> ■ tan (z/y), and r is the radius from the axis, r «* / y*+z . 
Cylindrical coordinates were chosen over Cartesian coordinates in order to avoid 
infinite body slopes in the cross-sectional planes. The body shape is generated 


from the coordinates (y^ A » z^) of cross sections at several axial stations x. , 


< 2 * ♦ • • * (see fig* 2) 


'ty -ij 


m 




Side view 


View from rear 


Figure 2. - Body coordinates 


Here a two-dimensional spline function is needed to describe the body, and the 
body shape may be a general three-dimens ion&l tody except .for the following 
restrictions: 

1. The body must have at least one plane of symmetry. 

2. The body radius r = f(x, <f>) must be single valued. 

3. The nose of the body must be located at x • 0 and at this position 
f(0, $) - 0 and |“ (0, <fr) ** 

Since the body has at least one plane of symmetry,, let this plane be the 
x - y plane and only the ha3.f body 0 < $ < it need be described. A simple doubly 
cubic spline (see Chapter VII of ref. 23) is used to generate the body shape 
from the input coordinates. The rectangular region 0 < x < x^, 0 $ <j> < ir is 
subdivided into a family of subrectangles x^_^ < x < x^, S <t>. S where 

0«x 1 < x 2 < . . . x^ and < .... ■ tt. Thus the simple doubly cubic 

spline for this two-dimensional mesh yields a function f^(x,$) which is: (1) an 

element of C 2 ^(R), where C r n (R) is the family of functions f A (x,<|>) on the rec- 
tangular region R whose n-th order partial derivatives , involving no more than 
r-th order differentiation with respect to a single variable, exist and are con- 
tinuous; (2) is a double cubic in each subrectangle, and (3) satisfies 
f <x ,* k ) * f^ for each point of the double mesh. 

To see how the one-dimensional cubic spline can be extended to the two- 
dimensional cubic spline, apply the one dimensional spline function of 

equation (5) to each cross section of the body at x^, x^> ...» x^. If the co- 
ordinates are chosen to be the same in each cross section, then the spline 

function becomes 


1^ 

6 + m j (x i^ 6 + 

4 1 ($**$,. 3) 
f J (x i ) - m J (x i ) 6 A 




n 1-l ( *i )A *1 

6 


A^>j 

(14) 


The quantities which change from one x-station to another are mj_^(x^), m^ (x^) , 
f J-l (x i)* ® nd f j ^ x i^ * Hence, could be applied to any general axial 

position to yield f^Cx,^) if the x-variation of m^U), m^ (x) , f^_^(x), and 
fj(x) were known. These variations may be determined by forming one-dimensional 
cubic splines of m^ (x) and f j (x) from the values of m^ (x^) and f ^ (x^) at 
x - x 1 , x 2 » ...» Xj^ for each value of j ■ 1, 2, ...» j m . 
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The one-dimensional cubic spline function f ^ (x) presents an immediate 

df (0) 

problem because the end condition *■ » at the nose cannot be handled. 

However, for the region 0 « x^ < x < x 2 , the axial variation of the radius 
r ■ f (x, <j>) for <J> ■ constant can be described by a general conic section for 
most bodies of interest. This general conic section may be written as 


[f(x,*)] 2 - 2R (<J>)x + [e 2 (*) - l]x 2 

V 


( 15 ) 


where R o (<f>) is the nose radius of curvature and e <4>) is the eccentricity of the 

general conic section (e ■ 0 for a circle, 0 < e < 1 for an ellipse, e ■ 1 for 

2 

a parabola, and e > 1 for a hyperbola). Therefore, if the function f (x,$) is 
used for the doubly cubic spline, in lieu of f(x,$), this function is well be~ 
haved at x « 0, because 


and 


Thus for the region 0 


which means that 


and also 


- ».<♦> 


3x 

x 1 < X < x 2 it is seen that 

- 2[e 2 (ij)) - 1) 
3x 


aVtaj,*) s 2 f 2 (x 2 ,^) 


3x 


3x 


(16) 

(17) 


(18) 


(19) 


3 3 4 f 2 (x^) 3 4 f 2 (x 2 ,$) 

2 2 " 2 9 

3x 3$ 3x 

These two equations are needed for the end conditions below. 

In view of the foregoing considerations, the doubly cubic Bpline 
2 

tation for f (x,<Ji) is 


( 20 ) 


represen- 
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U k - ♦)' 


(♦ - Vi y 


t^.t) - Vl («) -r^ + \< it > ~Ta V 

2 

9 Hr (fyr " ♦) 

+ !£ ( X i - m (x) -j-1 — — 

k 

2 

2 Afy, (<f> " \.i) 

♦ ' f k« - -6"> —^7- 


( 21 ) 


where the points 0 «■ <■ <p^ * . . • <()^ are chosen to be equally spaced around 

zn 

the circumference at each axial station (a^ “ ^ " $^-1 “ constant )* In 
equation (21) 

3 2 f 2 (x,$) 

\<x) - «— (22) 

* 3 <p 

2 

and m^(x) and f^(x) are determined from one-dimensional cubic spline fits to the 
data m^x^) and f^ x i^* respectively, for i ■ 1, 2, . i^. For x i _^ < x < 
these functions take the form 

(x 1 - x) 3 (x - x i 1 ) 3 

\ (1 ° • ~ir~ * B i.k e ax~ 


AX 2 (x - x) 
+ lm i-l,k " B i-l,k “6"! Ax^” 


Axf (x - x. .) 

- K, k - B i,k nr 1 —sr" 


(23) 


Although the input body coordinates are not generally equally spaced in the 

circumferential direction and generally differ from one body cross section to 

another, a one-dimensional spline fit to the data at each cross section can be 

2 

used to calculate new points f^( x i > <t>^) which are equally spaced in the cir- 
cumferential direction. 


and 


2 ( x .r " 

f k (x) “ A i _ 1<k + A 


<* - Vi r 


AX, 


i,k 6 Ax. 


2 Ax (x - x) 

+ _ i il i 

L i-l,k A i-l,k 6 1 Ax, 


2 Ax^ (x - X .) 

4. ff* _ a — — i.1 iZ A 

lt i,k A i,k 6 J Ax. 


( 24 ) 


where 


S 2 f 2 (x 1 ,«)> k ) 


L i»k = 

Sx 2 


sV <vV 

i»k " 

aq> 2 Sx 2 

III 

<3 

X i _ X i-1 


> (25) 


The end conditions at x ■ 0 follow from equations (19) and (20) as 


B 


l.k 


B 


2,k 


and 



In addition, it is assumed that at the end of the body 


(26) 


B i -l,k " B i ,k 
m m 


and 


V l ' k 


i ,k 
m’ 


(27) 


S 2 f 2 (x^.) 

These last two conditions imply that r — ~ and f (x,d>. ) have at most a 

k 

quadratic variation in x over the last subinterval [x.» _i * x^ ] . 

1 m 1 1 m J 

In the application of the cubic spline, it should be noted that the spline 
function passes through every input coordinate point. Therefore,, if an input 
point is incorrect, the spline function still passes through that point. On 
the other hand, the resulting derivatives from the spline function represent a 
smoothing of the actual derivatives (see pp. 42-44 of ref. 23). When accurate 
second derivatives are desired, a spline-fit of the first derivatives (called 
a spline-on-spline) has been found to be an effective tool. 
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Inviscid Surface Streamlines 


In order to apply the axisymmetrle analogue to three-dimensional boundary 
layers* it is necessary to trace out inviscid surface streamlines and calculate 
the corresponding "equivalent radius" or scale factor along each streamline, 

This scale factor (h) is the metric coefficient for the coordinate 6 normal to 
the streamline and on the surface of the body, It is a measure of the diver- 
gence of the streamlines as they wrap around the body. 

An orthogonal coordinate system will be used for the streamline geometry * 
Let the differential of arc length along a streamline be dS « h d£, the differ- 

B 

ential of arch length normal to a streamline (but on the body surface) be dq ■ 
h dg, and the differential of distance normal to the surface be dn. Both h 

6 

and h are metric coefficients or scale factors for the coordinates £ and g, 
respectively; however, only the scale factor h is needed for the application 
of the axisymmetric analogue. Since the coordinate n is the straight-line 
distance normal to the surface, its scale factor is unity » Along a streamline 
the coordinate g is constant, but it. varies from one streamline to another. 
Define the body geometry by Monge's form r ■ f(x,$) in a cylindrical co- 
ordinate system with the unit vectors e , e , and e in the x, r, and b direc- 

X t (p 

tions, respectively (see fig, 3). Note that the cross section at the right of 
figure 3 is viewed from the rear and that the coordinate system is right handed. 



Side View 


Rear View 


Figure 3. Body Geometry 
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Define 


F(x,r,<j>) 5 r - - C 

then the unit vector normal (outer) to the surface is 


Mi , * 1 8£ * 

; . n , iJxJklJj. ~.f 

n ]VF| 1/2 


( 28 ) 


where 


„ ,8f*2 f l 3f.2 

8 “ + 1 + If 8*1 


(29) 


Define the body angle f(x,<|>) by 


sin r * 


3f 

3x 


1/2 


<30) 


where -tt/2 < T < tt / 2 and also the body angle <S^(x,ij)) by 


sin 6 


1 11 

f 3J) 


[i Mf ff) 2 ] 1/2 


(31) 


where -tt/2 tt/2 (see fig. 4). From eqs. (30) and (31) it also follows 

that 

1 


cos r <« 


cos 6, g 

<P 


1/2 


and 


cos 6^ *• 


* [l+(j~) 2 ] 1/2 


Now the unit normal vector (eq. (28)) becomes 


e ■ - sin r $ + cos r (cos 6, e - sin 6 J ej 
n x <)> jr «J> 4> 


(32) 


As shown in fig. 4 the unit vector 


Sm ■ 008 6 * ®r ■ sln ®* (33 > 

is perpendicular to the curve of the body in a cross-sectional plane (but gener 

ally not normal to the body surface). Using eq. (33) in eq. (32), the unit 



* For an axisymmetric body, rr = 0 and hence 6 , = 0. 



normal vector can be written compactly as 


e - - sin f e + cos f e.,., (34) 

n x NN 

The unit vector tangent to the body curve in a cross-sectional plane (and also 
tangent to the body, see fig. 4) is 

i T i \ x e 'mi ‘ C0 ‘ 6 * S + sln S * e "r <35) 

and the unit vector tangent to the body surface and normal to e^ is 

•ll = ® n x ® T “ C0S F ®x + 8in r ®NN (36) 

Thus ® n » and e^, are a set of mutually perpendicular unit vectors with e^ 

and e^, tangent to the body surface. 



In order to determine the slope of an inviscid surface streamline, let 6 

s 

be a unit vector in the direction of the streamline with 8 the angle between « 

e g and e^ . Then 
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( 37 ) 


e = cos e e,, + sin 9 e_ 
s 11 T 

where 0 < 9 < tt and it is easily verified that e • & ■ 0. Define e as 

~ S XI g 

the unit vector tangent to the surface and normal to e g * i.e. 

e a = e x e = - sin 9 e,, + cos 9 a- (38) 

p s n 11 i 

Here also e , e , and e e are a set of mutually perpendicular unit vector with e 
SUP s 

and a tangent to the body surface. In terms of the unit vectors in cylindrical 

p 

coordinates, eqs. (33), (35), and (36^1 may be used in eqs. (37) and (38) to 

write e and e„ as 
s 8 

e = cos 9 cos F e + (sin 9 sin 6, + cos 9 cos <$ sin T)e 
s x 9 9 r 


and 


+ (sin 9 cos 6. - cos 6 sin <5, sin r)^ 

9 9 <1 


e„ « - sin 9 cos T e + (cos 9 sin 6. - sin T sin 9 cos 3,)e 
8 x 9 9 r 


(39) 


+ (cos 9 cos <S, + sin 9 sin 6. sin r)e. (40) 

9 9 9 

In order to trace out an inviscid surface streamline, it is first neces- 
sary to determine the angle 9 along the streamline. This angle can be ob- 
tained from the inviscid momentum equations (Euler’s equations) applied along 
the surface of the body. 

Euler’fi equation may be written as 

DV _ -Vp 
Dt ™ p 

where the operator is the time derivative along a streamline (for steady 

flows) and is related to the derivative D/DS as follows: 


JL a — — «» V D_ 

Dt “ Dt DS DS 

D 13 

where ~ » — — is the derivative along a streamline and again 
DS h o E 
s 

DS « dS ■ h d£ is the differential of arc length along a streamline. With 
6 

V d V e the acceleration vector becomes 



DV _ V DV ^ V DV 
Et " DS DS e s 


„2 De 

V s 


The pressure gradient in streamline coordinates is 

9p . 5E 4 * - - a -2- e + iL e 
p DS s h 86 6 an n 


Thus Euler's equations become 


V DV 


DS e s + 


„2 Da 

V S 


I$£e + 1 1£ e r l£e 

p DS s h dB ®6 an n 


The scalar product of e with eq- (43) yields the familiar equation 

S 

V DV 1 Dp 
DS “ ~ p DS 


where it is noted that 


S ^ ^ 

DS ' S ■ 0 


follows from e * e * 1, Next, the scalar product of e 0 with eq. (43) gives 

B S P 


De 
V s 


— — - • » 
DS 3 


1 dp 

Ph as 


To evaluate the left side of this equation, the scalar product of with the 
derivative of eq, (39) results in 


de • e. 


d8 + sin 1 do 


and the angle a is shown on fig, 4. Substitute eq. (46) into Eq. (45) to 


obtain 


M + 31,1 r M * ‘ “2 K ft (W 

This is the streamline equation in that it gives the angle 8 along a stream- 
line, However, before it can be applied, transformation operators are neces- 
sary to relate and ^ to derivatives in cylindrical coordinates. These 


* In obtaining eq. (46) from (39), note that de ■ e do and de, 

r p <j> 


d<p , 
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transformation operators are derived belowt 

Let 5 be the position vector of any point on the body surface relative 
to the origin of the cylindrical coordinate system (x « 0, r ■ 0) . Since the 
body normal coordinate n is zero on the surface of the body where r - f(x,<f>), 
the position vector can be written as 5 «* 5(x,<(>) and thus 


e dx + e, f dd> 

X <|> T 


(49) 


(see ref. 24). In terms of streamline coordinates on the body surface, the 
position vector can also be written as 


and thus 


$ - 5(5, $) 

d5 - || d5 + || dg 


where 


and 


■ e h d? + e a h dS 
s s P 


sl 


e h 
s s 


— - e h 

as n 

By equating the right side of eqs. (49) and (50), there results 


(50) 

(51) 

(52) 


e dx + e . f d<ji ■ e h d5 + e. h dft 
x T ss S 


(53) 


The scalar product of this equation with e^ yields 


dx » e • e h d5 + e. • e h dfi 

S X s $ X 


and since 


it follows that 


dx “ If d « + If d * 


3x a 

— MB t a n 

3 5 s x s 


(54) 
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Similarly, the scalar product of e^ with eq. (53) provides 

fdd> » e • e, h d( + e, • e h d6 

Y s <J> s p <!> 

and therefore 

t II - S • S. h (56) 

s 9 s 

and 

*H-vs h (57) 

Eqs. (54) to (57) may be used along with eqs. (39) and (40) to develop the 
transformation operators as 


1 8_ 1 1* £_ + I Hi- 

ll 38 “ h 3B 3x h 38 3* 


. 9 . V “( ) 

e * ‘ e p 57 + f a? 


^ (cos 0 cos 6^ + sin 0 sin 6^ sin r) ^ 
- sin 0 cos r 7— + 2 ~j 2 


(58) 


D 1 3 1 3x 3 , 134> 3 

DS “ h 35 “ h 35 Sx h 3~5 3«f> 
s s s 


« e 


^ V , 

e s 3x 



e 

s 


3_ 

3(j) 


(sin 0 cos <$. - cos 0 sin 6 sin r) 


■ cos 0 cos r — + 


3_ 

3<(> 


(59! 


In the equations above it is implied that all derivatives are evaluated on the 
body where r « f(x,(f>) and n - 0; thus 



I 

\ 


I 

k 


The transformation operators given by eqs. (58) and (59) can now be used to 
write the streamline equation (eq. (48)) in its final form as 


Si. . 

DS 


s 


p V 2 

8 


^ r 

, p .y 


- sin 6 cos r fj (|-) 


(cos 6 cos 6, + sin $ sin 6, sin r) . 

+ * 7 * It <?-> 


- sin r 


3cr 

cos 6 cos T -r— + 

0 K 


3<t> P, 


(sin 0 cos 6, - 


L 


cos 6 sin 6^ sin T) 
f " H 


(60) 


For a given surface pressure distribution, whether experimental or theoretical, 
this equation can be integrated numerically to determine 6 along a streamline 
emanating from the stagnation point. However, eq. (60) is indeterminate at the 
stagnation point, and this topic is covered below under the sub-heading Stag- 
nation Region Streamlines. The density p and speed V may be obtained from 
the pressure and entropy by use of isentropic relations for a perfect gas or 
equilibrium air. The geometric location of each integration step along the 
streamline is calculated by numerically integrating eqs. (54) and (56), which 
by the use of eq. (39) become 


and 


Dx 

ns 


cos 6 cos T 


M. 

DS 


sin 0 cos 6, - cos 0 sin 6, sin T 

1 

f 


( 61 ) 

(62) 


Since the surface streamlines and their corresponding orthogonal lines are 

curvilinear, the scale factors (metric coefficients) h and h play an im- 

s 

portant role because in the relation dS ■ h d£ the quantity d£ is an exact 
differential whereas dS in general is not. Likewise dg is an exact dif- 
ferential whereas the quantity dq » h dg is not. The importance of this state- 
ment lies in the fact that since d£ and dg are exact differentials, mixed 
partial derivatives are interchangeable, i.e. 


ail ail 

3538 ■ 3835 
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whereas it is noted that 


a 2 r , a 2 
asaq f aqas 

From eqs. (51) and (52) there results 


— ^ ^ = ^ / U L \ 

alas “ as as' as & 


ill » 1_ ll - i_ (h e ) 

asas 35 as as s s' 


By equating the right side of the two equations above, it follows that 


H (h V " a? (h s ®s^ 


which upon expanding becomes 


3e. 3h 3e 

111 i x v ! B — 1 e + h — — (63) 

as e s + h as as e s + b as 

The scalar product of e^ with the equation above yields 

Hi - h !i . a 
as a as s 

The equation above may be combined with eq. (46) and the relation dS ■ h^ dS 
to give 

511, ii + 8in r If. (64) 

ds as ae 

This equation cannot be used to calculate the scale factor h along a stream- 
line because ||- is unknown when only surface pressures are given. 

In order to obtain an equation for the scale factor h, rewrite eq. (48) as 

_ lE « ii + ain r 10 (48) 

oV 2 h 3 * 35 35 


* — : 3e ? 

As before, rr 


• e„ • 0 follows from e Q • e c ■ 1 



Differentiate this equation with respect to $ and substract the result from 
the derivative of eq. (64) with respect to £ to get 


1 dSl 

h DS 1 2 


i _ a_ 

h h 36 

s 



12 . 

36 


D(sin T) 3o 
DS h36 


Do_ 3(sin D 
DS h36 


where the equalities 


(65) 


2 2 2 2 
3 9 3 6 . 3 o 3 o 

35?6 “ 363? fin 3£36 " 363? 


were used to obtain eq. (65). This is the differential equation to be inte*- 
grate.d along a streamline to determine the scale factor h; however the terms 
on the right side of eq. (65) must first be cast into a usable form. 

By use of the transformation operators from eqs. (58) and (59) • it is 
found that the last two terms in eq. (65) may be expressed as 

D(sin r) 3o Do 3 (sin r) cos [ C0S jfc 

DS h36 " DS h38 " f 


3T. 3o 
3x a 4> 


3o ar 

3x 3 


( 66 ) 


and the first term on the right side of eq. (65) may be expanded to give 


1 3_ 

h h 36 
s 


JL.1E 


pV 2 h 86 


. 3h 
_1 s 

h h 36 


1 3p 
\ pV 2 h ^ 


1 3p 1 3o 1 3p 1 1 3V 

" h M p 2 v 2 h 36 " h 36 pV 4 h 36 


i_i a_ 
2 h 36 


PV 


112. 

h 36 


The four terms on the right side of eq. (67) are evaluated below. 

Take the scalar product of e with eq. (63) to get 

s 


3e s 

h — S. • 

3£ ' 


38 


(67) 


( 68 ) 
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Since e 3 . ^ . 0> then 


de 


de i 

ITr * e m _ s 

s rr * % 


and eq 8t (46) and ( 6 

y b ; combw •«. <«) to yl8l d 

ah 


- 

- h 

» 


ae i 

H + ,ln J'|f 




(69 


Comparing this equation with 8Q , ^ ^ (? ° 

(67 > becomes * * 16 ia eeen that the term -A. 

h h sg ln eq 


2 ah 

iTh rf * ~4- H 

8 pV h 33 


8 ^ pV"h 33 (71 

Since the inviscirf 

v 2 a ld fl W alon 3 the surface i* 4 

tSrm “h if in eq. (67) can be expressed a8 nt f°Pic, P - p( p ,s) and thi 


V. do 1 $ 

h 3s '(liphlt-M 2 A|E 

Thfi Bernoull i equation pV dV - - do 

in e S‘ <67) tc be written as 8 ^ SUrface H*** the term Ail 

h sg 


(72 


1 w v a v 


‘hw*-arff 

Iha fou tth and last term to t <73 

^ a PPiication o 'TJll: r 1 ^ <«, i. if fli*] 

i"S relation i, obtained: “°” 9per<it ° r o £ eq. <58 > the foU* 


f IT If I" fe . *ijLL iZ ] [ &] 

^ t h a 6 J [ ds I + s in e sin r Ab. A ar 

. J ax h al 

+ a i£ r 

f 3 $ sin e sin « cos f A 3T , , 

L * h 33 + (-cos 8 Sin 5 

<? 

+ sin 8 cos 6 Cos p n ,1 3fi * V e, 1 

* r> + .in 2 e co.Mi 

J . 2 


' 2 sine COS f (j, . i j l3j . , _ ? 

3 <P f 3x30 + (®g * A) 2 A» i_£ 

y a* 2 


8 


(74) 



In order to obtain eq, (74) , eq. (64) was used to substitute for r* 

n 3 p 

the term e^ * e^ follows from eq. (40) as 

e 0 • e. - cos 6 cos o , + sin 0 sin 6. sin r i 

P <P <p <p 

and eq. (58) was used to write 


^ gj ■ cos 0 sin <5^ - sin 6 sin r cos 6^ 


Using the results from eqa. (66) through (73), the equation for the scale 
factor h (eq. (65)) becomes 


2 

1 D_h 
h DS 2 


Ps p a V -l3 , , .I „ m 2v 

2 p 2 h 3$ (p/p s } (3 - M ) 

^ Q » 


, P s p s 1 3 13 . , . 

+ 2 p 2 h 38 I h 30 Cp/ V 

a m 


cos“ f coe 6 


3g _ 3o 
3<j> 3x 


Eq. (74) is used in the center term on the right side of eq. (75), and the 
transformation operators of eqs. (58) and (59) are used to transform partial 
derivatives from streamline coordinates to cylindrical coordinates. 

The geodesic curvature of a curve £ " constant on the body surface is 


and it is a measure of the amount that the streamlines converge or diverge. 

If ^ is positive h increases along the streamline, and. two neighboring 
streamlines move further apart, i.e. the streamlines diverge. Conversely, if 
the streamlines converge is negative. As indicated in ref. 6, the form of 
eq. (75) shows that h is not completely determinate. It may be multiplied 
by a constant or any function of 6. However, changes in h will also result 
in changes in 8. For the analysis here it is convenient to choose h and 6 
in such a way that they raduce to the radius r and the circumferential angle 
respectively, for the special cat* of an axisymmetric body at aero angle-of 


attack. Thus 6 will be dimensionless, and h will have the dimensions of a 
length. 

In summary, the inviscld surface streamlines and their corresponding scale 
factor h are calculated by numerically -integrating eqs, (60), (61), (62), and 
(75) for y, x, and h, respectively, along each streamline. The 

numerical integration scheme used here is the 4th order Venable step-size 
Runge-Kutta method. Initial conditions required to start the integration of 
each streamline are developed below in the section on Stagnation Region Stream- 
lines. See the Appendix for additional information on h 

It can be shown that eqs, ((>*) and (70) are the same as the Malnardi- 
Codazzi relations in ref. 24, and eq • (65) is equivalent to the Gauss charac- 
teristic equation in ref. 24 

Simplified Streamlines 

A simplified method of approximating the streamline direction on the body 
surface was used with good success by DeJarnette and Davis (ref. 14) in calcu- 
lating laminar heating rates over blunt-nosed cones at an angle of attack. In 
this method It is assumed that the direction of an inviscld surface streamline 
Is given by the resultant of the free-stream velocity vector minus its compo- 
nent normal to the surface*. This assumption was motivated by the fact that 
in the Newton-Busemann theory the initial direction of a particle, entering the 
shock layer is in the direction of these simplified streamlines. However, the 
Newton-Busemann theory has the particles following surface geodesics after 
entering the shock' layer, which differ from the simplified streamlines (see 
ref. 7), The simplified streamlines described here are called the method of 
steepest descent in ref, 7 and Newtonian streamlines in ref. 9. DeJarnette 
and Tei (ref. 15) found that the simplified streamlines could be used to cal- 
culate reasonably accurate- laminar heating rates for some cases even when the 
streamlines themselves were inaccurate, Simplified streamlines are developed 
here as an alternate method for calculating the inviscld surface streamlines 
and corresponding scale factor. 


The magnitude of the velocity, however, is calculated from the surface 
pressure and entropy, 
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Define as a unit vector in the direction of the free-stream velocity 

vector, then the simplified method gives the direction of a streamline by the 
equation 


e = •■ 

s i 


V - (V • e )e 
» ” n / n 

]y - (V • e )e 
‘ “ » n n 


e x (V x e ) 
n » n / 

1% * ( t x 


With V w in the body plane of symmetry (x - y plane) , it can be expressed as 
(see fig. 3) 

V w ■ cos a e^ - sin a (e f cos <j> - e^ sin p) (77) 

where the term in parentheses is the unit vector in the y-direction. Using eqs. 
(32) and (77) , the angle \p is given by 

cos i p 5 -V • e « cos o sin T + sin a cos r cos a (78) 

where 0 < < tt. It then follows that 


e x (V x e ) 
n 00 n 


V x e 
i » n 


sin i p 


and eq. (76) becomes 


V + cos \b e 

; b _!! L_n 

s sin <p 


To determine the streamline direction in terms of 8, equate the expres- 
sions for e given by eqs. (39) and (80) to get 
s 


sin 8 *» 


sin a sin o 
sin ip 


and also 


cor 8 - cos a " G0S It sin r . cos a cos r - sin a sin T cos q 
cos r sin ip sin \ p 

(0 < 8 < it) 


These equations along with eqs. (61) and (62) are used to compute the streamline 
geometry by the simplified method. In contrast to the previous section, which 
calculated 8 from a differential equation involving the surface pressure 


distribution, eq. (82) allows 9 to be determined at any position on the body 
without integrating a differential equation from the stagnation point to the 
point in question. 

It remains to develop an equation for the scale factor h corresponding 
to the simplified streamlines. Recall that eq. (64) was 


Dh se, r 3o 

DS SS + in f 39 


( 64 ) 


This equation can be used here to calculate h 
determined as follows. Differentiate eq. (78) to obtain 


36 


because the term rr can be 

dp 


3 T 

- sin ij» ■ (cos a cos T - sin a sin T cos c) gj 


3c 

- sin a cos r sin o 

0 p 


(83) 


and by using eqs. (81) and (82) this equation becomes 


- - CO. e |i + o: 


_ . . 3c 

cos r sin 6 


(84) 


Next, differentiate eq. (82) with respect to $ and substitute eq. (84) for 

i t 


36 


to -get 


If ■ fiT* <sln 8 la + “• 9 e0 ‘ r !f> ‘ sin r If 


Substitute eq. (85) into (64) to obtain 


_! Dh co s }J) 

h DS sin i}> 


1 3 r 

sin 0 h Tg + 008 6 cos r 


1 3c 1 

h 38 J • 


(85) 


( 86 ) 


1 3 


Then using the transformation operator of eq. (58) for rrr on the right side, 

n dp 

this first-order differential equation can be integrated aiong a simplified 
streamline to determine the corresponding scale factor- Note, however, that 

both eqs. (81) and (86) are indeterminate at the stagnation point, and the 

analysis in the section below must be used to supply the initial conditions for 
each streamline. Although this method of computing simplified streamlines 
should only be used on the windward side of the body, it was also used on the 
leeward side in reference 14 with reasonably good results for the corresponding 
heating rates. 


32 


Stagnation Region Streamlines 

The equations developed previously for the invlscid surface streamlines 
and scale factor h are singular at the stagnation point. In fact it will be 
shown below that this is a nodal point. Therefore, an analytic solution will 
be developed for a small region surrounding the stagnation point. This solu- 
tion will then provide Initial conditions to start the numerical integration of 
the differential equations for the streamline geometry and corresponding scale 
factor h given in the previous two sections. 

The actual location of the true stagnation point would require a numerical 
solution of the inviscid flow field over these three-dimensional bodies at an 
angle of attack. However, the true stagnation point is generally close to the 
Newtonian stagnation point for blunt-nosed bodies. Therefore, to be consistent 
with both the modified Newtonian pressure distribution and the simplified 
streamlines, the Newtonian stagnation point will be used in the analysis here. 
This point is determined by the position on the windward side of the body where 

V • e » - 1 
n 

This condition requires that r - - ct, and when ct > 0 symmetry dictates that 

6-0 and 6, » 0 at the stagnation point. 

For the region surrounding the stagnation point, it is advantageous to use 
coordinates 3^ and which are along the body surface and in the directions 
of e &-.d e^ , respectively. These unit vectors were defined by eqs. (35) 
and (36), and figure 5 illustrates these quantities. 



Figure 5.- Stagnation region coordinates 
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The coordinate S u is along the windward streamline (in the plane of symme- 
try), whereas is normal to S u but on the body surface. 

Consider again the position vector, relative to the nose of the body 

(x ■ 0, r ■ 0) as used previously. Then 5 can be considered a function of the 
streamline coordinates (5,8) or the surface coordinates (S^S^). It follows 
that 


^ - If « + f! ■ fs ^ 1 ds n + f!^ ds T 

■ h a d5 e s + hd& e^ - dS u e^ + dS T e^, (87) 

where it should be noted that dS^ and dS^ are not exact differentials. 

Take the scalar product of e^ with the equation above to get (using 

dS ■ h d£) 

s 

dS^ ■ dS e g • ejj + h dg eg • (88) 

Recall that *L j hence by using eq. (37) it follows from eq. (88) that 


DS 

dT 


11 


e^ ■ cos 8 


(89) 


In a similar manner it is found that 

ds t 

DS " ®s ‘ ®T " Sln 6 

which when combined with the previous equation yields 

DS_ 


(90) 


DS 


tan 8 


(91) 


11 


The inviscid surface velocity vector can be written as 

V — V e “V e + V_ e_ 
s 11 11 T 7C 


(92) 


Considering V n - V n (S u , S T ) and V T - V T (S n , Sj) , the stagnation point 
(S n “0 and S T » 0) requires that V n (0, 0) - V T (0, 0) » 0, and due to 
symmetry. 


3V n (s n* 0> 


98 . 


and 


v T (S n , 0 ) - 0 
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the stagnation point 

rherefore, *01 the region near the 
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stagnation point. Then eg. t»> — 

DS m 1 


T _ 


1 S T 


DS 


U 


B S 


11 


which may he Integrated to yield 

s 


C(fS) S 


B 


(94) 


(95) 


C(6) 1» really * £ut>cti0n ° J % * tnC * 

where the — - 

B8 t j_ u „ g la held conetant. 

the derivative Mn c0 „stant along a streem- 

tlnguiahea one streamline gndic.te. that the stagnation^ 

line, for convex bodies • * ttm region la spherical. » 

point is a nodal point. en stagnation point, for * 

the streamlines estate th . stagnation point, and the P 

streamlines do not emanate *•“»“» gained from eqs. (9» ond (»« 

the Streamlines at the stagnation p int (M) 

- fels-^lTicdr 1 


tane 8 - Ja " L* G(&) j 

fot c(B) finite and non-eero, the slope of a stre^ ^ ^ ^ whereas 
Hence, for q for B < 1 ot t&n 6 s 


motion point is tan 


for B « 1 (sphere) tan 6 g -*• which gives an infinite number of values as 
$ changes from one streamline to another. It will be shown later that by the 
use of modified Newtonian pressures B becomes the ratio of the two principal 
radii of curvature at the stagnation point' 

In order to evaluate C(3), consider the developed region around the stag- 
nation point, as shown in figure 6« 



Figure 6. Stagnation region 

On the circle of radius e, corresponding to the point where a streamline 
crosses this circle, the angle $ is assigned as shown in figure 6. Also, 
on this circle the coordinates of the streamline are S ■ e cos $ and 

B 

S T ■ e sin g. Now the streamline equation ■ C(3) is assumed to hold 

throughout this circular region. Applying this equation at the point where 
the streamline crosses the circle yields 

T5 

e cos 3 « C(3) (e sin 3) 

which gives 

C(3) - g ■ (97) 

(e sin S) 

Thus the parameter C(3) also depends on the value of e chosen. 
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The equation for the streamlines is now 



cos 8 
sin B 0 



( 98 ) 


Different values of 0 « 8 < it yield different streamlines. The slope of a 
streamline follows from eqs. (91) and (98) as 


tan 0 - 


5 

DS 


11 


ein B 0 , vl-B 

ITTosT (s t / e) 


(99) 


and for 8 the value of 6 on the e circle, where S_ 

£ * X 

equation above gives 


e sin 0, the 


tan 6 - (100) 

Hence for B 1 * 1, 0 g t 0; whereas for the sphere (B ■ 1) 8 g » 0 (radial 
streamlines). The windward streamline is 0 * 0 and the leeward is 0 ■ 180°. 

Next, consider the evaluation of the scale factor, h, in this region 
surrounding the stagnation point. Take the scalar product of e^ with 
equation (87) to get 


has - ds u 8 U • s 8 + as T Sj • 


Considering 6 « B(S . S T ), the previous equetion gives 


h H; ■ ' . % ■ cos 6 


Hence h » 


cos 6 

36 

3S m 


( 101 ) 


The derivative If’- is obtained by differentiating eq. (98), and the result is 
3S T 


80 d B cos 6 sin 0 

3S T S t (sin 2 0 + B cou 2 0) 

The term cos 6 may be obtained from equation (99) as 

B cos 8 


( 102 ) 


cos 0 ■ 


rt.2 2 a , , 2B 0 _ 28-2. 1/2 

(B cos 0 + sin & | s x| J 


(103) 
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Substituting eqs . (102) and (103) into eq, (101) yields the following equation 
for the scale factor 


h 


(sin 2 g + B cob 2 3)S t 


[B 2 cos 2 0 + sin 2B & 


(fT ) 2 ‘ 2B ) 1/2 sin 6 


(104) 


On the circle of radius e, where ■ e sin 3, the scale factor becomes 


, (sin 2 S + B cos 2 g)e 

* C.lA+B 2 co.V /2 


(105) 


For later use, is needed, and this quantity can be determined from the 


derivative of eq. (104). The result is 


fDh I _ 

. 2 2 12 
sin 8 + B cos 8 

Us u 

2 2 2 
. sin 8 + B cos 3 . 


(106) 


Again, it is worth noting that for a spherical region 
(106) give 


h 


e 


e 


and 



B 


1 


1 and eqs. (105) and 


which to first order are correct. 

As shown by Reshotko, ref. 25, the ratio of velocity gradients, 
evaluated from modified Newtonian theory as 



B, may be 


(107) 


where 


*T 


- 30 

008 r 8* 


(108) 


is the surface radius of curvature in the direction, and 


R 


11 


1 

cos r 


8T 

3x 


s 


(109) 


is the surface radius of curvature in the S,^ direction, both at the stagr 
nation point. 
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The location of each streamline on the e circle is determined by purely 
geometrical considerations. Suppose the stagnation region inside the e 
circle is represented by part of an ellipsoid with principal radii of curva- 
ture given by R^, and R^. Then for e/R^, << 1 the coordinates of a stream- 
line on the e circle are: 


x 

y 


e 


e 



+ e[cos 8 sin ot + — (sin 2 3 + B cos 2 ®)] 

+ e [ cos 3 cos a - (sin 2 3 + B cos 2 3)] 


( 110 ) 

( 111 ) 


s " e sin 3 
e 


(112) 


<J> e ■ tan 


-1 


(, A> 


(113) 


Note that the form of the above equations causes no difficulty at a * 0. 

In choosing a value for e, a compromise between accuracy and computational 
time must be made. The smaller the value of e the more accurate eqs. (110) 
through (113) become. On the other hand, it was found that the step si2e used 
to start the numerical integration of the streamline equations outside the e 
circle must be less than about e/10. After testing several values of e, a 
value in the range 0.01 < e/R T < 0.1 was found to be sufficiently small to 
make eqs. (110) through (113) reasonably accurate, yet large enough to keep 
the integration step size for the streamline equations from becoming prohi- 
bitively small. 


Surface Pressure Distribution 

As mentioned previously, the calculation of the invlscid surface stream- 
lines and heating rates are dependent on the surface pressure distribution. 

Two options are considered here for the pressure distribution, depending on 
whether pressure data from another source is available or not. First, if 
experimental or theoretical pressures (ratioed to stagnation pressure) are 
known, tabulated values around the circumference at several axial stations 
are used to generate a two-dlmensicnal spline function for the pressure ratios 
and derivatives at any position on the body. This two-dimensional spline 
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function is similar to that used to describe the body geometry except for the 

axial variation. For blunt-nosed bodies, !§■■*■» at x ■ 0 but is 

354 3/x 

finite there. Therefore, the pressure spline function Is developed to vary 
with /sT in the axial direction. For a given circumferential angle <J> the 
derivative for both x ■ 0 and the end of the body is determined by the 

derivative from Lagranglan interpolation. 

When surface pressures are not available from some other source, the 
second option is to use modified Newtonian pressures, l.e. 

2- - (1 - ~)cos 2 * + — (114) 

Pg Pg Pg 

where the angle i|> is given by eq. (78). Since derivatives of the pressure are 
needed to calculate the streamline geometry and scale factor, .modified Newtonian 
pressure ratios are computed for each body coordinate used to generate the body 
geometry, and then the same pressure spline function described above is used to 
spline fit this data. In the "shadowed" region (cos ^ < 0) , p » p w is used. 

Gas Properties at Edge of Boundary Layer 

In the calculation of heating rates, the local pressure (p g ) , density (p g ) , 
enthalpy (h ) , velocity (U ) , speed of sound (a ) , and coefficient of viscosity 
(y ) are needed at the edge of the boundary layer. The first-order boundary- 
layer approximations allow the pressure at the edge of the boundary layer to be 
the same as the corresponding surface pressure (p * p) , which was discussed 
in the previous section. Also, the boundary layer is assumed to remain suf- 
ficiently thin so that the entropy at the edge of the boundary layer is constant 
and equal to that value aft of a normal shock wave. This section describes how 
to calculate the other flow-field properties from the pressure and stagnation, 
properties. 

Equilibrium air properties . - Many elaborate computer programs are avail- 
able for calculating equilibrium air properties. However, these more sophisti- 
cated approaches require considerable computational time and storage locations. 
In order to keep thr computations relatively simple, the correlation formulas 
of Cohen (ref. 22) re used here to calculate the equilibrium air properties. 
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2 


\ 


These formulas were shown to be quite accurate in ref. 22 for pressures in the 
range 10 ^ atm. <_.p < 10 atm. and an enthalpy range from 128.7 BTU/lb (corre- 
sponding to a temperature of 540°F) to 16,930 BTU/lb (corresponding to flight 
at approximately 29,000 ft/sec). 

Assuming an isentropic expansion from the stagnation point, eq. (65) in 

ref. 17 gives the enthalpy (h ) as 

6 



h tl - 2.119 x 10 8 ft 2 /sec 2 


(115) 


This approximate equation for h is more restrictive than the limitations 

given previously. Beckwith and Cohen (ref. 17) found that eq. (115) gave 

results generally within 3 percent of highly accurate computerized properties 

except when p > 5 atm with H <• 4,500 BTU/lb. 
s s 

With h given by eq. (115), the density can be computed from the equa- 
tlon (see ref. 22) 


5 „ * 965 

7.344 x 10“ 3 (y*jy) 


* 

1 - 1.0477 1 - | 

fh \.6123 
e 

m 

\h 1 


(116) 


and then the coefficient of viscosity is 


2.0144 x 10~ 10 (y^y) 


.992 




1.0213 


1 - 


.33291 


i 


(117) 


The velocity is obtained from the adiabatic energy equation as 

». ■ - v ]1/2 


( 118 ) 



! 
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Although the local speed of sound a @ Is not given in ref. 22 , it can be 
readily Calculated as follows. Since the entropy at the edge of the boundary 
layer Is assumed to be constant , 


d E_ 


<iv*t 


The inverse derivative dp Q /dp can be. obtained from the derivative of eq. 
(116), using the isentropic relation dh * dp/p . The result is 



Finally, the Mach number follows from 


M 


e 


u 

e 


( 120 ) 


( 121 ) 


Perfect gas properties. - For a perfect gas with y ■ 1.4, the usual 
isentropic relations, as given in ref. 27, may be used to obtain the following 
equations : 


h 


e 





( 122 ) 


P 


e 



(123) 


U e - [2(H g - h e )] i/2 


(124) 


T 

e 


h 

m e 

6006 


(125) 


The coefficient of viscosity is computed from Sutherland's law, which gives 


2,27 x 10~ 8 T 3/2 
e 

y e “ ~ T + 198.6 
e 


(126) 
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and the speed of sound is simply 


a g ■ (1. 4)p/p e (12' 

Stagnat i o n pro perties* - In the analysis here the stagnation streamline 
is assumed to pass through the normal part of the bow shock wave. It is then 
necessary to calculate both equilibrium air and perfect gas stagnation-proper' 
ties from known values of p , p , T , and V . 

First, consider equilibrium air. The conservation equations for mass, 
momentum, and energy across a normal shock wave are: 


V » " P 2 V 2 


+ P. K " P 2 


+ P. V * v 2 


H 

s 


K + 



09 


2 


h 2 + 



(121 


(13< 


These equations must be solved iteratively for equilibrium air, and the pro- 
cedure used here is as follows: 

1. Assume » 0. 

2. Calculate p 2 from eq, (129). 

3. Calculate h 2 from eq. (130). 

4. Calculate p 2 from eq. (116). 

5. Compute a new value of v 2 from eq. (128); then go back to step 2 
and start the process again. 

6. Repeat steps 2 through 5 until the new value of v 2 is sufficiently 
close to its previous value to give convergence. For an accuracy of 
10 8 , six iterations will generally suffice. 

Using p 2 , h 2 , and H g , eq. (115) can be rearranged to give 


(it h - E 

'2117 ‘ 

.035 

h„ 

A 

.3877 

/ H , 

f 8 

.3877 

1 ,1 

v.035 

\ 3.345 x 10 8 

>7\ 




. l h 2 1 

1 

I 


( 13 : 

then eq. (116) can be used to calculate p . 
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whsr© the factor 778 has been added to make q have the dimensions of 

2 w 

BTU/ft -sec, This equation can be applied independently to any inviscid 
surface streamline. However, it becomes indeterminate at the stagnation 
point since both and h go to zero there. This indetermancy can be re- 
solved by the use of the results obtained in the section on Stagnation Region 
Streamlines. 

From eq. (104), the scale factor h in the stagnation region is 

2 2 
(sin S + B cos B) S m 


h « — 


G sin g 


where 


G £ [B 2 cos 2 g + sin 2B & (S T /e) 2_2B ] 1/2 


Replacing V by the velocity at the edge of the boundary layer U , eqs. (91) 
and (93) may be combined to give 

V. G S T <»V»V. 


sin 9 " . B„ . 1-B 

sin B (S T /e) 


and also 


“ a T 
sin 9 


G PS T 

sin B g (S T /e) 1_B 


Thep ' equations are accurate to the first order in the stagnation region, and 
to be consistent the pressure ratio to first order of accuracy is 

" 1 (139 

Eqs. (135) through (139) can now be used to evaluate the integral in eq. (134) 


, (sin 2 S + B cos 2 g ) 2 e 2 " 2B S 2+2B 


■ 2"R 

s sin - ' g (2 + 2B) 


which can be used to evaluate the indeterminate part of eq. (134) as follows: 
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(141) 


Using this result, the limiting form of eq. (134) becomes 
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v.-Mv® »*‘ ,M v/ 


u e >. 


3 V. 


3S„ 


(142) 


Consistent with Lees' approximations for a cold wall, H g “ H g - h^, 

(o y ) * (p y ) , and thus eq, (142) is essentially the same as Reshotko's 

K e e s w w s n 

equation (ref. 25). 

For the present analysis eq. (142) is modified so that it will reduce to 
two-dimensional and axisymmetric stagnation-point heat-transfer rates that are 
compatible with experimental results. In ref. 17 the following equations were 
found to compare reasonably well with experimental data: 
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Consistent with the analysis of ref. 26, the factors 0.577 and 0.768 in eqs. 
(143) and (144), respectively, are replaced by 


I 



(145) 


for the general three-dimensional stagnation point. Using the axisymmetric 
analogue, the pressure-gradient parameter 8 for a relatively cool end 
isothermal wall is determined from eq. (33) in ref. 17 as 
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where t g - h g /H g . For e sphere g g ■ 0,5, whereas if « 1 for a cylinder. 
Using eqs. (146) and (135) through (139), it is found that jf g for a general 
three-dimensional stagnation point is 


and 


6 . 

1 

" B+l 

for 

B 

> 1 


B 

“ B+l 

for 

0 

< B < 1 


1 

i (147) 
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However, Beckwith (ref. 5) found that g » 1 - | for .5 < B < 1, which gives 

® * 

the correct limiting values for a sphere (B ■ 1, g ■ ~) and a cylinder (B - 0, 
— s * __ • 

g g » 1) but disagrees with eq. (147) above. Eq. (147) gives g g « ~ for a 

sphere (B ■ i) , but for the cylinder (B - 0) it gives jf ■ 0. This result is 

s 

not surprising for the axisymmetric analogue because the streamlines are assumed 
to originate from a single stagnation point, and therefore all the streamlines 
on the cylinder are forced to emanate from this stagnation point rather than 
from the stagnation line along a generator of the cylinder. As will be shown 
below, the pressure-gradient parameter g has only a small influence on the 

o 

magnitude of the term given by eq. (145). 

From eqs. (52) and (59) in ref. 26, the following equation is obtained 
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Using this result in eq. (145), the modified form of eq. (142) consistent with 
eqs. (143) and (144) is 
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For a relatively cool wall the equation of state gives 


P w “ I 2 " (154) 

R T 

w 

Then for an isothermal wall and u * u (T ) , it follows that 

w w w 


By substituting eq, (141) 
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into (146) 
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, the following equation results. 


(155) 


The factor 



8 t 

e 


l 

[A] 

(B+l) 

s 

D 

f V 


V 

90 1 

AJ 

_DS 

U. J 

S 



-* 


T T 

S/^w,s in eq. (153) is obtained from eq. (61) in ref. 26 as 
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As suggested by Cohen (page 33 of ref. 26), the first factor on the right side 
of eq. (157) was used to replace the factor 1.033^ in eq, (61) of ref, 26. 
Since the wall is assumed co be relatively cool and isothermal, 



(158) 


The final form of the laminar heating-rate ratio is obtained by substitut- 
ing eqs. (146), (155), and (156) into eq. (153), which gives 
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where the factor ; /; Is obtained from eq. (157), B is determined from 

w w, s r 


eqs. (107) to (109), and 
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may be obtained from the modified Newtonian 


value given by eq. (152). 

Laminar heating rates are obtained by applying eq. (159) along an inviscid 
surface streamline from the stagnation region to the end of the body or the 
beginning of the transition region, whichever comes first, Keating rates along 
each streamline are computed independently of the other streamlines. The 
integral in the denominator of eq. (159) can be evaluated by quadrature for 
each integration step along a streamline. 


Transition Region Heating Rates 


The location of the transition region on a body is a subject open to much 
speculation and debate. It is not proposed that the present method will predict 
the location of the transition region. Instead, the beginning and end of the 
transition region may be specified by one of the three options listed below. 

1. Geometric, location, or 

2. A specified value of the integrated unit Reynolds number along a 


surface inviscid streamline, i.e. Rn 1 • 


P u 

H e e 


DS , or 


3. A specified value of the momentum thickness Reynolds number, 


p U e 

Rn „ _® — ® — SL | a i on g the inviscid surface streamline, 

m y e 

Once the beginning and and of transition is determined by one of the .three 
methods above, the heating rate in the transition region is calculated as a 
weighted average of the local laminar and turbulent heating rates. Thus the 
heating rate is written as 
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However, when transition is specified by the integrated unit Reynolds 
number or momentum thickenss Reynolds number, the corresponding geometric 
position along a streamline for the end of transition may lie off the body sur- 
face. In this situtation the value of S„ cannot bs determined. To circum- 

tre 

vent this situation, the "weighting" function for these two options will be 
based on the value of the integrated unit Reynolds number (Rn^) or momentum 

Therefore eq. (164) becomes 


thickness Reynolds number (Rn ) 

TU 


w^ * 1 - exp 


-.412 


4,74 (R^ - Rn lttrl > 
(R &l,tre " ^ljtri^ 
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(165) 


when transition is based on the Integrated unit Reynolds number, or 


w. 


1 - exp < -.412 


4,74 (R ^m “ Sutri* 
^ Rn m,tre " Rn m,tri^ 


-1 2 


(166) 


for transition based momentum thickness Reynolds number, where the subscript 
"tri" refers to the beginning of transition and "tre" refers to the end of 
transition. 


Turbulent Heat-Transfer Rates 

Cooke and Hall (ref. 6) have shown that the axisymmetric analogue is 
applicable to turbulent as well as laminar boundary layers. There are many 
methods that could be used to calculate turbulent heating rates on axisymmetric 
bodies. The approach used here is to apply the axisymmetric analogue to a 
modified form of the integral method of Reshotko and Tucker (ref. 20) to obtain 
the turbulent momentum thickness. Then the momentum thickness is used to cal- 
culate the local skin friction coefficient from the correlation formula of 
Spalding and Chi (ref. 19), and finally the Karman form of Reynolds Analogy 
(see ref. 21) is used to calculate the turbulent heating rate corresponding to 
the skin friction coefficient. 

Momentum thickness. - Reshotko and Tucker (ref. 20) developed an Integral 
method to calculate the momentum thickness from the numerical solution of two 
first-order differential equations for 0 and H^ along a streamline for a 
perfect gas. This method is modified here to apply to a perfect gas or 
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equilibrium air, and the skin friction coefficient used in the method is that 
value computed simultaneously from the Spalding-Chi approach described below. 

The basic method used by ref. 20 can be applied to both perfect and equi- 
librium gases if the ratio a /a is replaced by A /H , For the perfect 

gas these two expressions are one and the same; however, they are different for 

1 DM 

an equiliblrum gas. As a result, the quantity rr- must be replaced by 

H s Bu e e 

k u DS~ * T ^ en t * le inte 8 ral * orm the momentum equation becomes 
e e 


De m 2. ®m D ^e ®m Dh C f 

ST* (2 + H - M e ) — 
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where is the local skin friction coefficient. Following ref. 20, the 

compressible form factor is given by 


where 
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The incompressible form factor is calculated from the equation 
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(169) 


[2.6 (H t - 1) - -43 (H 1 + l) 2 ] - .03^ 0 l ± - 1) y 


(170) 


where due account has been made for a perfect gas or equilibrium air. 

Eqs. (167) and (170) are integrated numerically along an inviscld surface 
streamline, starting at the beginning of the transition region. The initial 
values of 6^ and H are the laminar values given by ref. 17, and C f is 
supplied by the Spalding-Chi approach described below. 

Skin friction coefficient . - The turbulent skin-friction coefficient is 
determined by the correlation formula of Spalding and Chi (ref. 19) baaed on 
momentum thickness Reynolds number but modified to account for perfect or 
equilibrium gases. Eq, (52) of ref. 19 is modified to read 
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is the turbulent recovery factor. 
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With the definitions 

(173) 


U* - (2/C f F c ) 1/2 (174) 

the compressible .form of eq. (28) in ref. 19 becomes 

F rs R { » f «S> 2 + (KE)' 1 ^ U - (2 /KUq) ] exp <OT*) 

t (2/KuJ) + 1 - | «£) 2 - 15 «U+) 3 

» I— (KU + )^ — (KU + )^1 (175) 

40 ^ G ; 180 V *V J 

whera K ** 0*4 and E M 12* In tha application of ©q* (1/5)* and &r© 

computed from eqs, (171) and (173)* respectively* with 0 obtained from the 
Reahotko-Tucker method described above. The Uq Is calculated from ®Q* (175) 
by an iterative scheme such as the Newton-Rhepeon method. Using U Q , eq. (174) 
gives the incompressible skin friction coefficient since 


Flnallyj 0^ 
In ref. 


follows from eq. (176) once F £ 
19. F is defined by 


has been evaluated. 


( 176 ) 



where p and U are the density and velocity across the turbulent boundary 
layer. The boundary layer approximations give p « p , and for the perfect gas 
the equation of state yields 


where h^ is the enthalpy inside the boundary layer, which is assumed to 
follow Crocco's relation 
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Upon substituting eqs. (173) and ,(179) into (177), the resulting expression can 
be Integrated for the perfect gas to yield 
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Tor equilibrium air eq. (180) is not valid because eq. (173) is based on 
a perfect gas. The density-enthalpy relation for equilibrium air is obtained 
from eq. (116) as 


£ 


1 - 1.0477 [1 - (h e /h E ) ,6l? ’ 3 } 


p e 1 - 1.0477 [1 - (h e /h E )‘ 6123 (h bl /h e )* 6123 ] 


(181) 


Then using eq, (179) in. (181), the integral in eq. (177) is evaluated numeri- 
cally to obtain F c for equilibrium air. 

Turbulent Heating-Rate Expression. - The Kerman form of the Reynolds 
analogy factor is given by ref, 21 as 
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e Pr 0.725 was used. As in ref. 21, it is assumed chat the compressible 
cyno s ana ogy factor is the same as the incompressible factor, i.e. 


Then, since 
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the turbulent heating rate follows from eqs, (182) -(184) as 
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Both C f and C f>1 are obtained from th. Spalding-Chi approach described above. 


COMPUTATIONAL METHOD 

The method developed herein was programmed on the IBM 360/75 digital 
computer at th. North Carolina State University. This program is also com- 
pat . with the CDC 6600 computer at th. Langley Research Center, NASA. A 

Z tnZT* # * the C °" PUt “ * PPe " S in P “‘ 11 - — 

To run a typical case, the following data are needed as input parameters 

T >- v »> «. v 

2. specify a perfect gas ( Y - 1.4) or equilibrium air. 

3. body geometry data as Cartesian coordinate, of points around the 
circumference at several axial stations. 

4. specify surfac. pressure distribution by one of the following three 
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a) p/p from some other source for points around the circumference 
at several axial stations along the body, or 

b) modified Newtonian pressure distribution, or 

c) modified Newtonian pressure distribution, but streamlines 
computed by the simplified method of ref. 14. 

5. specify the beginning and end of the transition region by one of the 
following three methods: 

a) geometric location, which is specified by the Cartesian coordi- 
nates of points around the circumference of the body for both 
the beginning and end of transition, or 

b) values of Rn^ at the beginning and end of transition, or 

c) values of Rn^ at the beginning and end of transition. 

If transitional and turbulent heating rates are not required, lami- 
nar heating rates alone may be calculated. 

6. The value of 8 for each inviscld surface streamline to be computed, 
where 0 a 0 < 180°, 8 “ 0 is the streamline in the windward plane 
and 8 ** 180° is the streamline in the leeward plane. 

All the input data used to calculate the results in this part of the 
report are given in Part II. The computer program calculates heating rates and 
other pertient data along each inviscld surface streamline Independently of the 
other streamlines. 


RESULTS 

To illustrate the present method, results are presented for blunted cones 
at angles of attack, a blunted 70° slab delta wing, and the HL-10 lifting body. 

Blunted Cones 

Blunted 9 B half-angle cone . - Figure 7 shows three inviscid surface 
streamlines (8 K 5 s , 20 s , and 90°) calculated for a blunted 9° half-angle cone 
at a ■ 10° and » 18 using a perfect gas (y * 1.4). Each of the three 
streamlines was calculated by three different methods: (1) method of charac- 
teristics (ref, 28), (2) present method using modified Newtonian pressures, 
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and (3) simplified streamlines by the method of reference 14. Although the 
geometry of this body could be represented by simple analytical expressions, 
the coordinates of 19 points around the circumference of the half body at 20 
axial stations were used to generate the bray geometry by the doubly cubic 
spline function. Figure 7 shows that the & * 5 1 and 20 s streamlines calculated 
by the present method using modified Newtonian pressures are in good agreement 
with the "exact" streamlines calculated by the method of characteristics. The 
@ * 90* streamline calculated by the present method does not agree well with 
the "exact" result, but this is to be expected since this streamline goes into 
the "shadowed" region (<t> ± 154*) on the conical afterbody where the pressure 
is set equal to free-stream static pressure. The simplified streamlines are 
significantly lower than those calculated by the other two methods. 

Blunted 15; half-angle cone. - In reference 29 Cleary gives tabulated 
experimental laminar heating-rates on a spherically blunted 15° half-angj.e cone 
at M w ■ 10.6, In order to compare results from the present method with some of 
this data, streamlines and laminar heating rates were calculated on this con- 
figuration for a » 20°, R q - 0.375", p w « 2,6614 lb/ft 2 , T w « 89,971°R, and 
V w ■ 4928,1 ft/sec, (M ■ 10,6). Although a perfect gas with y • 1,4 was used 

in the calculations, gas imperfections in the wind tunnel stagnation properties 
5 2 

of p “ 1,73 x 10 lb/ft ana T « 20Q0“R were taken into account to determine 
s s 

the free-stream conditions. In addition, it was determined that a value of 

X, S h /H ■ 0.251 corresponds to T /T * 0, 270, 

*w w s r w a 

The body geometry was specified by che coordinates of 20 points around 
the circumference of the half body at 19 axial stations. Streamline patterns, 
calculated by the present method using modified Newtonian pressures and by the 
simplified method of ref. 14, are shown in figure 8 for 8*1®, 10®, 15® and 
45®. This figure shows that the streamlines move from the windward plane and 
rapidly approach the leeward plane, and the streamlines calculated by modified 
Newtonian pressures wrap around the surface at a steeper rate than the simpli- 
fied streamlines. The "equivalent radius" or scale factor increases along a 
streamline on the windward side, where the streamlines are diverging, but 
decreases along the leeward side due to converging streamlines. As mentioned 
previously, the streamlines are continued into the "shadowed" region of the 
body C<t> 2 137.4°), although the pressure is assumed to be free-stream static 
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pressure in that region- It is possible that the flow separates from the sur- 
face somewhere In this region, and therefore the- calculated streamlines and 
heating rates are questionable there . 

Laminar heating-rate ratios (q /q„ ) , calculated by the present method 

using both the streamlines computed from modified Newtonian pressures and the 
simplified streamlines, are compared with Cleary's experimental data* in fig- 
ures 9, 10, and 11, Figure 9 shows the heating-rate ratio along the windward 
plane, whereas figures 10 and 11 show the circumferential distribution of 
heating-rate ratios at axial stations of x ■ 3 >56" and 9.36", respectively 
(corresponding to x /L • 0,207 and 0,466 in ref, 29). Very good agreement 
with the experimental data is obtained for the heating-rate ratios using 
modified Newtonian pressures and reasonably good agreement using simplified 
streamlines for this case. For this relatively simple body shape, the heating- 
rate ratios calculated using the simplified streamlines are close to those 
using modified Newtonian pressures although the streamlines themselves differ 
by a large amount (see fig, 8). It is significant to note that the ratio of 
the angle-of-attack to the cone half-angle is 1.33 for this case, and thus the 
present method is not limited to small angles of attack. 

Blunted 70® Slab Delta Wing 

This configuration is a 70® swept delta wing with a cylindrical leading 
edge which is tangent to a flat slab on the upper and lower surfaces. The 
blunt nose is a spherical c'ap, and the radius of both the spherical cap and 
cylindrical leading edge is one foot. Although this configuration could also 
be specified analytically, coordinates of 13 points around the semi-periphery 
at 19 axial stations (0 £ x s 1.95 ft) were used to generate the body shape. 

Figure 12 illustrates the streamline patterns calculated by the present 
method (using modified Newtonian pressures and a perfect gas) on this delta 
wing at ct ■ 10® and ■ 8 (p^ « 10^ lb/ft^, T^ » 416®R, V w ■ 8000 ft/sec, and 
5 ® 0.4), It can be seen that the 8 « 10® streamline is converging towards 

the center line on the flat slab which causes the scale factor to start 

The experimental heating rates in ref. 29 are ratioed to the calculated 

2 

q„ ■ 35,94 BTU/ft -sec in this report. 

W| 8 
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decreasing. On the other hand, the 6 * 60“ streamline runs nearly parallel to 
the leading edge beyond the nose, and at x « 2 ft the scale factor is still In- 
creasing which indicates the streamlines in this region are diverging. The 
streamline patterns in figure 12 are qualitatively similar to the experimental 
oil-flow patterns in reference 30 for this angle of attack. 

Laminar, transitional, and. turbulent heating-rate ratios along the wind- 
ward streamline (<j> « 0) are shown in figure 13- Transition was arbitrarily 

.. ~ 

chosen to begin at Rn • 10 and end at Rn - 2 x 10 . These values were used 

m m 

to illustrate the capability of the present method to calculate transitional 
and turbulent heating races, and they should not be construed to represent-the 
actual transition region on this body. As shown in figure 13, the transition 
region corresponding to these values of Rn^ lies in the narrow band 
0.268 ft < x « 0,318 ft, This figure also shows chat the heating-rate ratio 
increases sharply in the transition region, and then decreases in the fully 
turbulent region- For x i L ft the windward streamline is on the flat slab 
where the heating-rate ratio decreases very slowly. 

HL-LQ Lifting Body 

The geometry of the HL-10 lifting body without fins is illustrated by plan 
and side views in figure 14 and cross sections in figure 13. This body is an 
example of a shape whose geometry is difficult to describe analytically. In the 
present method the body geometry was generated from the coordinates of 20 points 
around the semi-periphery at 20 cxial stations. 

Figure 16 depicts laminar, transitional, and turbulent heating-rate ratios 
calculated along the center line of the lower surface for o *» 20®, M w ■ 10, 
p - 10 3 lb/ft 2 , T » 416 S R, V » 10 5 ft/sec, and 5. ■ 0.1. Equilibrium air 
and simplified streamlines were used in these calculations, and transition was 
arbitrarily chosen to begin at x *» 33,07" and end at x » 54.24". Again, the 
heating-rate ratio increases sharply in the transition region, and then begins 
to decrease near the end of that region- 

Difficulties were encountered for this case when modified Newtonian pres- 
sures were used to calculate the inviscid surface streamlines. Along the wind- 
ward streamline the scale factor went to zero, which indicates merging or 
crossing of streamlines and invalidates the heating-rate calculations. This 
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difficulty can be traced to the fact chat away from the nose the lower surface 
of the HL’-10 body is flat in the z-direetion, and modified Newtonian pressures 
are constant across a flat surface.. As shown in the Appendix , this feature of 
modified Newtonian pressures generally produces questionable results « 

DISCUSSION 

In the calculation of inviscid surface streamlines and scale factors , 
first and second derivatives of both Che body geometry and pressure are needed. 
The doubly cubic spline function was found to represent the geometry of the 
bodies considered herein satisfactorily- However, its accuracy is affected by 
the number and location of body positions used, and the spline function does 
not smooth out data points (although it does smooth derivatives; . When pressure 
data from some other source are used in the present method, they should be 
smoothed before using them so that the pressure spline function will be 
accurate. 

Inviscid surface streamlines calculated from modified Newtonian pressures 
compared well with those from the method of characteristics for a 9* half-angle 
cone at a » 10° and M w * 18, except for the "shadowed" region. In the 
"shadowed" region on the leeward side, the pressure is assumed to be constant 
at p M , and thus the streamlines calculated in this region are Inaccurate and 
follow geodesics of the surface (see ref. 7), In addition, the flow may sepa- 
rate from the surface somewhere in the "shadowed" region, which makes the 
results from the present method questionable there. 

The present method was found to predict laminar heating rates very well on 
a blunted 15° half-angle cone at a *> 20* and M w * 10.6 using modified Newtonian 
pressures. Additional comparisons of the present theory with experimental data 
on other body shapes and for transitional and turbulent heating should be made 
to assess the accuracy of the theory more thoroughly. 

Bodies with flat segments, like the bLunted 70 s slab delta wing and the 
HL-.10 lifting body, present difficulties when inviscid surface streamlines are 
calculated using modified Newtonian pressures or simplified streamlines.. On a 
flat segment, both the modified Newtonian pressure and the angle c are constant, 
and consequently eq. (48) shows that the streamlines are straight (6 • constant), 
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and the direction of each streamline is determined by its direction upon enter- 
ing that flat segment. Since each streamline may have a different direction at 
the beginning of the flat segment, it is possible for the streamlines calculated 

from modified Newtonian pressures to cross over one another*. In addition, eq. 

2 2 

(75) shows that D h/DS » 0 along these calculated streamlines which means that 
the scale factor h is llnearily increasing or decreasing along the streamline, 
depending on the value of Dh/rs. at the beginning of the flat segment. In par- 
ticular, if Dh/DS < 0 it is easy to see how h calculated by this approximate 
method may go to zero (as encountered on the HL-10 body) and even become nega- 
tive. The scale factor along the windward streamline is discussed more thor- 
oughly in the Appendix. On a flat segment the simplified streamlines are all 
parallel and the scale factor h along each streamline is constant and equal 
to the value the streamline had upon entering the flat segment. The difficul- 
ties encountered on flat segments with these two approximate methods for calcu- 
lating lnviscid surface streamlines can be circumvented by calculating the 
streamlines from an accurate pressure distribution. 

In the present method the boundary layer is assumed to be sufficiently thin 
that the flow may be considered to be isentropic at the edge of the boundary 
layer. This assumption may be invalid away from the stagnation point on some 
body shapes due to the boundary layer "swallowing" the entropy layer. When 
this occurs, the entropy at the edge of the boundary layer could be signifi- 
cantly less than the entropy aft of a normal shock wave. 

The computational time required to compute the examples given herein is 
highly dependent on the number of streamlines specified in the input data. 
Typical cases with six streamlines calculated from modified Newtonian pressures 
required less than 4 minutes on the IBM 360/75 computer for a perfect gas. When 
simplified streamlines were used, the computational time was reduced about 50%. 
On the other hand, the computational time was increased about 10% when equilib- 
rium air was used in place of a perfect gas. 


Of course, it is physically impossible for one surface streamline to cross 
another. 
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CONCLUDING REMARKS 


A method is developed for calculating mviscid surface streamlines and 

laminar, transitional, and turbulent heating rates on general blunt-nosed three- 

dimensional bodies at angles of attack in hypersonic flows- Relatively simple 

techniques are employed to keep the computational storage and run time down to 

small values (less chan 4 minutes for a typical case on the IBM 36C/75 computer) , 

Streamlines calculated from the modified Newtonian pressure distribution 

were found to compared favorably with those from the method of characteristics 

on the windward side of a blunted 9* half-angle cone at a. » 10' and M * 18, 

» 

Laminar heating rates calculated on a blunted 15° half-angle cone at oi » 20 s 
and M w * 10.6 compared very well with experimental data,, 

Streamlines and laminar, transitional, and turbulent heating rates were 
calculated on a blunted 7Q V slab delta wing at a » 10 s and M^ « 8, and on the 
HL-10 lifting body at u * 20' and M x • 10. These bodies have some flat seg:- 
ments, and since the modified Newtonian pressure is constant over a fiat sur- 
face, streamline patterns and corresponding heating rates calculated from 
modified Newtonian pressures are questionable. This difficulty can be circum- 
vented by using an accurate surface pressure distribution in lieu of modified 
Newtonian pressures. 

The relatively small amount of numerical computations required coupled 
with reasonably good accuracy, makes the present method attractive for engi- 
neering applications. Additional comparsions with experimental data on other 
body shapes and for transitional and turbulent heating are needed to assess 
the accuracy of this method more thoroughly. 

Mechanical and Aerospace Engineering Department, 

North Carolina State University, 

Raleigh, North Carolina, Augus.t .25,. JL971. 


APPENDIX 


Scale Factor Along Windward Streamline 

Along the windward streamline <J> ■ 0, 6 - 0, and symmetry requirements 
dictate that 6 - 0, 6 -0, dV/H * 0, and ^ \ • Therefore, for this 

streamline 
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is the body radius of curvature in a cross-sectional plane (x ■ constant) at 
iji « 0. When eqs, (Al) and (A2) are substituted into eq. (64), there results 
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This equation can be integrated to give the scale factor along the windward 
streamline only when the velocity-gradient term on the right side is known « 
For an axisymmetric body (o = <j>) , ■ f and eq« (A5) reduces to eq« (B-xl) 
in reference 11. 

When only the surface pressure distribution is known, eq. (75) must be 
used in lieu of eq. (A5) to calculate h, and along the windward streamline 
it reduces to 
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Figure 15, - Frontal cross sections of HL-10 lifting body. 




re 16 , - Heating-rate ratio along windwar- 



